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Abstract  -  For  many  oceanic  studies,  it  is 
required  to  know  the  distribution  of  visible  solar 
radiation  (EPAR)  in  the  upper  water  column. 
One  way  to  reach  this  is  by  remote  sensing.  This 
includes  two  components:  First,  EPAR  at  surface 
is  calculated  based  on  atmosphere  properties 
along  with  the  position  of  the  Sun.  Second,  the 
vertical  attenuation  of  EPAR  (KPAR)  is  derived 
from  products  of  ocean-color  remote  sensing. 
Currently,  KPAR  is  estimated  based  on 
chlorophyll  concentration  ([C))  from  ocean 
color.  This  kind  of  approach  works  well  for 
waters  where  all  optical  properties  can  be 
adequately  described  by  values  of  [CJ,  but  will 
result  in  large  uncertainties  for  coastal  waters 
where  [C]  alone  cannot  accurately  describe  the 
optical  properties.  In  this  paper,  we  present  an 
innovative  model  that  describes  KPAR  as  a 
function  of  water’s  inherent  optical  properties 
(IOP). 

I.  INTRODUCTION 

Solar  radiation  in  the  visible  domain  (EPAR(350 
-  700  nm),  measured  by  downwelling  irradiance  in 
this  text)  encompasses  the  wavelengths  shorter 
than  700  nm.  The  pioneer  study  of  Zaneveld  et  al. 
[1]  and  subsequent  studies  [2-4]  have  demonstrated 
that  the  vertical  penetration  of  EPAR  plays  an 
important  role  in  heat  transfer  of  the  upper  water 
column.  Epar  at  surface  can  now  be  adequately 
estimated  from  satellite  measurements  of 
atmosphere  properties.  It  requires  information  of 
water’s  optical  properties  to  determine  the  vertical 
attenuation  of  EPAR  (KPAR)  with  depth.  Historical 
measurements  have  shown  that  KPAR  not  only 


changes  horizontally  with  constituents  in  the  water 
[1,5],  but  also  changes  with  depth  for  any  water  [6, 

7]- 

To  represent  the  steeper  than  exponential 
reduction  of  EPAR  with  depth,  multiple  exponential 
terms  [6,  7]  were  usually  adopted,  with  an 
attenuation  coefficient  (or  attenuation  depth) 
assigned  for  each  term.  These  attenuation 
coefficients  are  kept  vertically  constant,  but 
horizontally  vary  with  Jerlov  [5]  water  types. 
Recently,  simple  and  explicit  models  have  been 
developed  to  incorporate  satellite-derived 
chlorophyll  concentrations  ([C])  into  the 
description  of  the  attenuation  of  EPAR.  When  [C] 
values  are  provided  via  satellite  observations  of 
ocean  color  [8,  9],  the  partition  factors  and 
attenuation  coefficients  of  the  terms  could  be 
calculated  [4]. 

Such  kind  of  approach  works  for  Case-1 
waters  -  where  all  optical  properties  are 
determined  by  [C]  alone  (with  solar  zenith  angle 
explicitly  or  implicitly  included)  [10,  11].  For  non- 
Case- 1  waters,  uncertainties  arise  due  to  that  it  is 
not  a  constant  relationship  between  [C]  and  optical 
properties.  To  avoid  such  limitations  associated 
with  [C] -based  models,  another  approach  is  to 
describe  the  vertical  transmittance  of  EPAR  using 
water’s  optical  properties  [12,  13].  Following  this 
strategy,  and  because  that  water’s  absorption  (a) 
and  backscattering  coefficients  (/ih)  can  be 
adequately  derived  from  ocean-color  remote 
sensing  [14-16],  we  developed  a  model  to  describe 
the  vertical  transmittance  of  EPAR  using  values  of  a 
and  bb. 
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II.  Hydrolight  SIMULATED  KPAR(z) 

As  in  earlier  studies  [17,  18],  we  used 
HydroLight  [19]  to  get  the  necessary  data  sets: 
Kpar(z),  and  a  and  Z>b.  Unlike  the  simulations  in 
Ohlmann  and  Siegel  [18]  where  water’s  IOPs  were 
determined  by  [C]  only,  IOPs  in  our  simulations 
were  simulated  with  varying  [C]  and  independently 
varying  CDOM  and  suspended  sediments,  as 
described  in  Lee  et  al.  [14]  and  IOCCG-OCAG 
[20].  Later,  KPAR(z)  is  modeled  as  a  function  of  a, 
bb  and  z.  Numerous  descriptions  can  be  found 
regarding  simulations  by  HydroLight.  [19,  21-23]. 
The  following  summarizes  the  input  settings 
carried  out  in  this  study. 

The  downwelling  irradiance  at  sea  surface 
from  the  Sun  and  sky  is  simulated  by  the  spectral 
model  of  Gregg  and  Carder  [24],  a  and  bb  values 
at  440  nm  varied  from  0.02  to  1.9  m"1  and  0.002  to 
0.115  m'1,  respectively,  and  kept  vertically  constant. 
The  wavelengths  are  in  a  range  of  350  -  700  nm 
with  a  10-nm  spectral  resolution.  Five  depths 
(excluding  0  m)  were  selected  for  each  HydroLight 
run,  with  depths  spread  within  and  beyond  the 
euphotic  zone  [25].  No  bottom  reflectance  and 
inelastic  scatterings  (such  as  Raman  scattering)  are 
included  in  this  study. 


III.  MODELING  OF  KPAR(z) 


With  Epar(z)  simulated  by  HydroLight 
Kpar(z)  is  calculated 


K 


par(z)-~  ln| 
z 


A7W0-)A 


- PAR 


(z) 


(1) 


Figure  1  presents  a  few  examples  of  KPAR(z). 
Clearly,  KPAR(z)  differs  significantly  for  varying 
water  properties.  Also,  consistent  with  earlier 
measurements,  subsurface  KPAR(z)  changes  a  lot 
even  for  vertically  homogeneous  waters.  This 
change  is  due  to  that  water  molecules  absorb 
strongly  in  the  longer  wavelengths  (large 
absorption  coefficients).  After  photons  pass 
through  the  subsurface  layer  (say  3  meters),  the 
absorption  is  happened  in  the  shorter  wavelengths, 
where  absorption  coefficients  are  generally  smaller, 
especially  for  oceanic  waters. 

For  each  vertical  variation  of  KPAR(z),  it  is 
found  that  this  vertical  change  could  be  modeled  as, 


Kpar(IOP,z)  =  K1(IOP)  + 


K2(IOP) 


,0.5 


(2) 


(1  +  z)u 

Here  K\  is  for  the  asymptotic  value  at  greater 
depths,  with  K2  more  important  to  the  subsurface 
Kpar  value.  IOP  here  represents  different 
combinations  of  absorption  and  backscattering 


coefficients.  The  dotted  lines  in  Figure  1  show 
Eq.2  modeled  Kpar(z)  for  those  examples.  Figure  2 
presents  the  result  of  HydroLight  KPAR(z)  versus 
Eq.2-modeled  KPAR(z),  with  the  Sun  at  30°  from 
zenith.  Apparently  the  modeled  KPAR(z)  matches 
the  Hydrolight  KPAR(z)  very  well  (the  average  error 
is  2.2%,  with  maximum  error  of  is  6.4%).  Such 
results  clearly  demonstrate  that  Eq.2  is  adequate  to 
describe  the  vertical  change  of  KPAR(z). 
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Figure  1.  Examples  of  KPAR(z)  for  different  water 
properties.  The  numbers  in  the  box  are  values 
of  a(490)  (left)  and  bh( 490)  (right).  Symbol 
represents  KPAR(z)  from  HydroLight 
simulations,  while  dotted  lines  are  models 
from  Eq.2. 


Figure  2.  KPAR(z)  from  model  (Eq.2)  compared 
with  Kpar(z)  from  HydroLight  (30°  solar 
zenith  angle),  indicating  that  Kpar(z)  can  be 
well  described  by  Eq.2  with  two  parameters. 


To  apply  IOP  distributions  obtained  from 
satellite  observation  of  water  color,  how  Kl  2  vary 
with  IOP  needs  to  be  known.  For  the  KPAR(z)  data 
with  the  Sun  at  30°  from  zenith,  it  is  found  that  Ki 
and  K2  could  be  well  modeled  with  IOPs  at  one 
wavelength  Zaneveld  et  al.  [12]  Barnard  et  al.  [13]: 


absorption  and  backscattering  coefficients  at  490 
nm,  i.e., 

K\ C rop )  =  Zo  +  X\  («(490))0'5  +  j2  bh (490) ,  (3a) 

K2(IOP)  =  g()  +  a(490)  +  ^,(490) .  (3b) 

2o,  i,2  and  go, 1,2  are  model  coefficients. 

Since  KPAR(z)  also  varies  with  solar  altitude, 
we  carried  out  HydroLight  simulations  with  the 
Sun  at  10°  and  60°  from  zenith  in  order  to  include 
solar  zenith  angle  into  the  model.  From  these 
simulations,  we  got 
K\(IOP,Ga)  = 

iZo  +  Zl («(490))°-5  +  j2  **(490)] (1  +  a0  sin(0a ))  ’ 

(4a) 

K2(IOP,0a)  = 

2  a  .  (4b) 

[fo+fia(490)  +  f2^(490)](a1  +a2cos(Oa)) 

Here  6a  is  the  solar  zenith  angle  above  the  surface. 
Now  we  got  a  model  that  can  describe  the  vertical 
distribution  of  EPAR  for  different  IOPs,  depth,  and 
sun  angle  as 

T(IOP,z,0a)  =  EpAR (z)  =  e~KpjR(IOP’z'da)z  .  (5) 
epar(®  ) 

In  this  model,  there  are  nine  model  coefficients: 
2o,  i,2,  go, i,2,  and  £*0,1,2-  To  derive  their  values,  T 
values  from  Eq.5  were  fit  against  T  values  from 
HydroLight  simulations  with  the  model  coefficients 
derived  by  least-square  curve  fitting  [4,  18].  Values 
of  derived  2o, 1,2,  go, 1,2,  and  £*b,i,2  are  provided  in 
Table  1.  Figure  3  presents  Eq.5  modeled 
T(IOP,z,Qa )  versus  TUOP.z.  6„)  determined  from 
HydroLight  simulations.  For  those  T  values 
(limiting  to  the  range  of  -0.001  to  0.8),  bigger 
errors  happened  at  T  <  0.003,  where  the  effects  of 
E? AR  on  heat  transfer  and  photosynthesis  in  the 
water  column  are  small.  For  T>  0.003,  the  average 
error  is  -9%.  These  results  indicate  that  the  simple 
optical-property-based  model  (Eq.5)  is  adequate 
for  describing  the  vertical  profile  of  EPAR(z)  for 
different  waters. 

IV.  SUMMARY 

In  this  study,  an  innovative  model  is  developed 
for  describing  the  vertical  transmittance  of  visible 
solar  radiation  (EPAR(350  -  700  nm),  measured  by 
downwelling  irradiance  in  this  study)  in  the  upper 
layer  of  the  oceans.  Different  from  the  traditional 
approaches,  one  exponential  term  is  used  for  the 
vertical  distribution  of  E?AR.  Its  attenuation 
coefficient  (KPAR(z)),  however,  is  modeled  as  a 
function  of  depth  instead  of  vertically  constant. 
Kpar(z)  is  also  modeled  as  a  function  of  solar 
zenith  angle  and  water’s  optical  properties  (a(490) 


and  bb(490))  with  data  from  HydroLight 
simulations.  With  the  availability  of  a(490)  and 
bb(490)  images  obtained  from  satellite  remote 
sensing,  this  KPAR(z)  model  can  be  adequately 
incorporated  into  physical  oceanography  models  to 
study  the  effects  of  visible  solar  radiation  on 
surface  heating  [26].  Also,  it  provides  easy  and 
reliable  tool  to  predict  the  light  level  at  desired 
depths,  needed  to  plan  the  C14  incubation  for  in  situ 
measurements  of  primary  production  [13]. 
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T  =  Epar(z)/Epar(0-) 


Figure  3.  T(z)  from  model  (Eq.5)  compared  with 
T(z)  from  HydroLight  for  three  sun  angles. 


TABLE  1.  MODEL  COEFFICIENTS  FOR  KPAR(z) 


parameters 

values 

Xo,  Zl?  %2 

-0.057;  0.482;  4.221 

?o;  <gi;  ?2 

0.183;  0.702; -2.567 

oc0;  oq;  0C2 

0.090;  1.465;  -0.667 
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